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ABSTRACT 


The design, construction and calibration of swmmeie 
instrumentation for measuring temperature fluctuations 
within a few centimeters of the sea surface, surface wave 
variations, and wind velocity are discussed. Small bead- 
in-glass thermistors mounted at 2.5 and 5 centimeters 
above the water surface were used to reasure temperature 
variations under varying wind conditions and spectra] 
densities computed on a digital computer. Spectra. co- 


spectra and coherence squares from 0 to 5 Hz are presented. 
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l, introduction 

The energy exchanges between the atmosphere and the 
oceans are the driving force for a thertodynaric and mech- 
anical system which encompasses the fluid portions of the 
earth. The atmosphere derives energy from the sea in the 
form of sensible and latent heat and in turn, provides 
mechanical energy for oceanic circulations. Small-scale 
processes, such as convection and local wind fields, are 
also direct results of these energy exchanges, and the 
ability to accurately measure energy fluxes at any level 
over the sea would permit one to predict the growth and 
decay of such processes, 

There are generally four forms of air-sea energy 
exchange that one may consider: 

Vane LeCeromasnel tc veadrar rom 

2. Convection and conduction of sensible heat 


S42 SCCOMmMectwon and dirruscton ol latent eames (eine 
form of water vapor 


4, Vertical exchange of nechanical energy. 
Iiiietnwestneating the flux of any of these forms of energy, 
ome DUeEwimcasure directly certain quantities. “Por heat 
flux, the appropriate equations are: E. = eC Tw and E, = 
LW, whe re i is the sensible heat flux. Ey is the latent 
heat flux,e anda, are the densities of air and water re- 
spectively, os is the specific heat of air at constant 
pressure, L is the latent heat of evaporation (or conden- 
sation). T is the temperature and w is the vertical cor- 


ponent of the wind velocity. Then, to measure the heat 


fluxes, one must be able to measure temperature, vapor 


pressure, and wind velocity at any desired level. Since 


waves on the water surface interact with the windeivemd 


(which carries the scalar temperature 


field), the sur- 


face wave field may also have an effect on heat flux, 


As a preliminary step in the determination of heat 


flux, methods of measuring variations of temperature, 


wind velocity, and water surface displacement on the 


smallest possible time and space scales must be devised 


and tested. This study was made to investigate some of 


the possible approaches to these measurements close to 


the air-sea interface and had as its objectives: 


1, The design of a suitable sensor system for 


measuring the temperature fluctuations within 


a few centireters of the sea surface, the wave 


field at the same point, and the downstream 


wind velocity at a suitable 
surface, 

24, The @ontinvieite recordimnever 
urements in a form suitable 
ysis. 

3, The analysis of the records 
appropriate spectra and the 


between the air-temperature 


level above the 


these three rmeac- 


for digital anal- 


to “determine reme 
relation, if any, 


VArFLactLOns and “the 


water-surface height under varying wind con- 


Gverens. 
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i: Desten of the Measurdhent Systen 
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A system to measure srail temperature variations 
requires a temperature sensor that is durable, inexpen- 
Sive, and that has a sufficiently small response tire to 
record fluctuations at frequencies of about 10 Hz. A 
small bead thermistor was used and although thermistor 
time constants are generally starge compared to, say, a 
platinum resistance thermometer, it appeared that they 
were sufficiently small for this investigation. A Vic- 
tory Engineering Corporation bead-in-glass thermistor, 
model 32A49, was used. The manufacturer's specifica- 
POMS EOr this Unit are: a nominal diameter of 0.0137, 

a time constant (7) of 1.0 seconds in stili air, a pow- 
Pematssipation constant of 0.10 mw/°C, and a continuous 
Operational range of -50°C to 300°C. The zero-power 
Peemstance at 25°C is 25002 a 29%, and the resistance- 
temperature ratios are 5.6 for Ks at O°C/R, at 5 On Gearand 
14.6 for R, at 25°C/R, at 125°C, However, over the range 
of temperatures encountered in this study (about 5°C), the 
resistance-terperature curve is essentially linear (Figure 
i . 

Two of these thermistors were mounted on a guided 
styrofoam float, at about 2.5 and 5 cm above the water 
surface, as shown in Figure 2. A variety of float designs 
were tested for response te surface waves. The one used 
in the field was a hemisphere of about 2.5 cm diameter. 


In wave tank tests, movies showed the float to follow the 


waves quite well, although, of course, it tended to act 
as a filter for most waves smaller than the diameter of 
the float. 

To measure the water-surface fluctuations, a resist- 
ance wave gauge was mounted as close as possible to the 
thermistors. This*gauge consists of a pair of 9201) ai 
ameter stainless-steel wires, about 25 cm in length and 
positioned some 3 cm apart (Figure 1). A 5 volt, 20 Hz 
power source was used to avoid electrolysis, and the sig- 
nal output from the bridge was fed through a full-wave 
rectifier. 

Wind velocity was measured with a Casella sensi- 
tive cup anemometer, mounted approximately one-half meter 
above the water. The anemometer has a digital counter to 
record wind velocity and was used to determine mean wind 
speeds over time intervals of 10 to 30 seconds, 

The wave gauge and therrmistors were installed as 
one leg of separate Wheatstone bridges and the outputs 
recorded on Mosley /l100B strip recorders. This recorder 
has individual DC amplifiers for each channel with maxi- 
mum amplification providing a 5-mv per 10-inch scale. 

The recorders were operated at 0.5 inches-per-second”@mart 
speed for all measurements. The bridge and rectifier cir- 
cultry for each instrument 1s” shown ier eurece se oncdee 
and a block diagram of the basic data-collection system 


is given in Figure 5. 


Lies 


The platforr for #ield=tounting them sensors is 
shown installed for measurements in Figure 6, The trans- 
misstom lines ror Signal output are shielded, two-core 


telemetry cables with low-loss characteristics. 


3; Instrument Calibration and Observation Techniques 
The thermistors were calibrated by varying the tem- 
perature of a water bath and plotting the resulting re- 
eieeeance Change Tor EG€ach Lemperature increment. the re- 
sults are presented in Figure 8. Since this curve was 
essentially the sane for each individual thermistor, it 
was assumed that minor differences in slope could be ig- 
nored. Before field measurements, the bridges were bal- 
anced in the laboratory at a temperature close to that 
expected in the field. Each recorder channel was cali- 
brated to a scale of 1°C per inch and for equal sensi- 
tivity in each channel. The entire system was then test- 
ed with a known temperature change with the telemetry 
cables installed. The response time of the system was 
estimated by running the recorder at its maximum chart 
speed (2 inches per second) and recording a known termper- 
ature change. The tire for 63% of the change (one time 
constant) was then reasured directly from the record. 
Since the temperature changes were induced by adding 
water of a different temperature to the bath, the actual 
temperature change was not instantaneous, and the esti- 
mate of 80 msec is perhaps too large. It is, however, a 


good indication that a flat response to 10 Hz is a con- 


eee 


servative estimate. This response will be somewhat less 
in moving air, and a response of 5 Hz was considered to 
be the upper limit for a flat response for this study. 

The wave gauge was calibrated in a wave tank to find 
the required immersion for maximum signal amplitude. The 
smaller the immersion, the greater the signal amplitude. 
Above a wave amplitude of about 5 cm, the output became 
slightly distorted due to the non-linearity of the gauge 
system. Again, the range of teasurements in this study 
was sufficiently small that the response could be con- 
sidered linear. The output of the full-wave rectifier 
was compared with a common generated signal of O to 10 Hz 
and showed an undistorted output over this range. 

The Casella anemonveter was not calibrated. The manu- 
facturer's conversion curves for counter readings to veloc- 
ities were assumed to be accurate, 

The field observations were made during the months 
of April and May, 1968 at Roberts Lake, Seaside, Califor- 
nia. This small hake is approximately one-half mile in- 
land from Monterey Bay and is about 450 yards long and 
200 yards wide. Between the bay and the lake are a series 
of sand dunes about 25 to 30 feet high which end about 150 
yards from the lake. There is no other obstruction to air 
flow in this area to seaward and the winds are essentially 
those present over Monterey Bay. The available umnrestrict- 
ed fetch at the measurement site is between 75 and 400 


yards, depending on the wind direction. Most of the data 


was taken with a fetch of from 100 to 300 yards. The 
platform was placed in about three feet of water, 25 to 
30 feet from shore, and the signals were trans~itted to 
the recording equipment with the low-loss cables. 
Measurements were made on six occasions, the first 
three with only a single thermistor rounted at 2.5 cm. 
The lengths of the recordings were between 3 and 5 min- 
utes, with mean measured winds from 2.4 to 8 meters per 
second. The significant wave heights were between 2 
and 4 cm, with periods of from 1.2 to |1.8 seconds and 
lengths of from 10 to 15 cm. One representative record 
of the single-probe measurerents and one of the dual- 
probe measurements were analyzed. A portion of a compu- 


ter drawn record of dual-nrobe data is shown in Figure 8. 


4. Analysis 

In previous air-teuperature measurements (Pond et 
al, 1966), a substantial portion of the spectral "energy" 
was concentrated in the region below 10 Hz. Since 5 Hz 
was considered to be near the upper limit of the flat 
response of the thermistor, the analysis of the data was 
restricted to the range of frequencies from 0 to 5 Hz. 
As the data was recorded in a graphic form, the conver- 
sion to digital form required a sanual interpretation at 
intervals of 0.05" and 0.025°C, This process certainly 
Maieroduced some nolse into the record (in addition to 
that generated by the electronics in the system), but this 


was unavoidable. The ensuing analysis showed this noise 


bs 


to be negligible. 

Since the temperature fluctuations at the air-sea 
interface are felt to be essentially randor and station- 
ary, it iS appropriate to analyze the data by computing 
power spectra. This method cf analysis is the basis of 
an existing FORTRAN IV program in the library of the 
N.P.G.S. computer library, PROGRAM BLACKY (Table I) 

This program computes smoothed power spectral estimates 
fror the Fourier transform of the covariance for a given 
frequency bandwidth. 

In analyzing any continuous process usirg discrete 
points, the probler of aliasing is encountered, For this 
study, it was assumed that iittle energy existed above 5 
Hz. PROGRAM BLACKY also proviced a filter by running a 
weighted average for threé consecutive points of the 
record. This would reduce the possibility of folding and 
also rerove most of the high frequency noise. 

In order to confirm the r@sults of PROGRAM BLACKY, 
and to obtain a graphic display of the spectra, the data 
was also run in PROGRAM BIMED O2T. an autocovariance and 
power spectrum program of the Health Science Facility, 
UCLA. The records were aisc analyzed for means, density 


aie cumulative distri out vous. 


2. _ Resuits_and Conc iusions 
The data was processed for szcothed estimates of 
spectral density trem 0 to 5 Hz@et en tater va see 


cycles. The record was fiitered by means of a Hanning 


ee 


lidiewindow of the form: 
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MN 
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DS CoE ote i 


where P is the weighting function, m the number of lags 
fieommtrme wecord, and tr is am incremental index of mae The 
weights for this analysis were 0.25, 90.50 for the three 
consecutive points. The first and last points were com- 
puted with weights of 0.5 a Ole x4 and 0.95 Xnep t 0.5 
Xn This low-pass filter effectively removed the high 
Frequency noise of the record, Since little energy was 
felt to exist at high frequencies, the loss of this por- 
tion of the signal was considered to be insignificant. 

Ihe smoothed spectral estimates for the temperature 
variations are presented in Figures 9 and 10, and the 
corresponding wave amplitude spectra in Figure ll. A co-= 
spectrum for the dual-probe temperatures is presented in 
Figure 12. The values fron which these plots were taken 
are given in Tables I! and III. The spectra are presented 
as log, ¥(£) versus logip(£), Mabislauees 4 (icy) eilygMeGnghnics (oe 


(°c)*sec, The spectrum is defined such that 
oo : 
{¥ (at = 
oO 


where 6 is the temperature variation (Batchelor, 1959). 
The wave amplitude spectrur, @(f), is defined similarly, 
with units of cm“sec. Since Taylor's hypothesis of "fro- 
zen turbulence” (Taylor, 1938), is not valid for a region 
of shear flow, the relation 


k = U 


pant 


-— 


where k is the wave number, { the frequency, and U the 
mean wind velocity, was not used and the spectra are pre- 
sented in terms of frequency rather than wave nucber. 

The total "energy" was estimated in two ways in PRO- 
GRAM BIMED 02T; (1) by computing the variance, and (2) by 
integrating the power spectruy with the trapezoidal rule. 
The ratio of <1)/(2) gives the fraction of the total ener- 
gy in the computed spectra. The following table gives the 


individual values and the resulting ratios. 


<< ome eer er ee ee ee ce ceo © ces ce ce em em mm SS a ee ee ee eee eee ee se ee 


cE y/OCE, Check sum Energy by Drrn. Ratio 
ee OT BS Dire a pil emreenenemeee nate ae 
Temp. 2.5 cm 1.1057844 1.1058750 ~-9.06x107? 1.00 
(Single-probe 
Temp. 2.5 cm 0.2300338 0,2300585 2.47x10°? 0.99 
(Dual-probe) 
Temp. 5 cm 0,3806309 0.3806635 3.27x107? 0,99 
(Duai-probe ) 

4 


Wave Amplitude “225749149 2 75751486 2.347310" 5 woe 
(Single-probe ) 


Wave Amplitude 0.2900739 2.2900738 1.00x107’ 1.00 
(Dual-probe) 


The high percentage of energy in the corputed spectra indi- 
cates that there is little energy above 5 Hz and there was 
negligible energy folding from higher frequencies. 
Confidence limits for the spectral estimates were 
calculated for an 80% interval. For the analyzed record 
length of 120 seconds, a sampling interval of 0.1 seconds, 


and a correlation lag value of 2C0, the 80% lizits are; 


ee 


0.65 ¥(f) = V(t) = 1.91V7(£). The covputatrion of these 
limits is presented in Appendix A, 

Coupling, between the surface waves and the tem- 
perature fluctuations at fixed heights above the undu- 
lating surface, does not appear to exist, even at the 


higher wind speeds. The coherence square function is 


2 an) = 
LE | 3k aera 
XY Gif) CG (£) 


where er ae is the co-spectrel density function and CkOs) 
and GC £) are the individual spectral density functions 
(Bendat and Piersol, 1966). The coherence square for the 
temperature and wave spectra (Table II) shows virtually no 
correlation between these quantities. Since the tempera- 
ture measurements were taken at a height that remained 
constant with respect to the water surface, this lack of 
correlation suggests that the temperature structure at 
this level is not affected by the passing wave surface. 
In contrast, spectral densities for temperature measure- 
ments taken at a fixed height above mean sea level have 
shown periodic maxima at the saizre frequencies as wave 
amplitude spectral maxima (Makova, 1963). A comparison 
of measurements taken at a fixed height above the surface 
and simultaneously at a fixed height above mean sea level 
would be useful in further studies of this problem. 

The one-dimensional spectra of temperature fluctua- 
tions plotted on a iog-log scale (Figures 9 and 10) show 


a ve slope for frequencies above 0.2 Hz. This slope has 


Lie 


also been found for temperature seasurerents taken with a 
platinum resistance thermometer at heights of one meter 
and greater (Pond et al, 1966). Since the inequality, 

kz 27 4.3. does, not hold near 0.2 Ha, this result» adndicanrec 


Ehatt a me lation of tChewmrerm 
iy > / 
Vif) = G Aenea Ss 
te 


S > 


(1) 


(Batcheicr, 1959) is valid fer anctielously large sca! 


(Here, C, is a conversion constant to relate the spectrum 


1 
to frequency rather than wave nucber, eg is the sceawar 
dissipation, e is the energy dissipation, and A is a uni- 
versal constant.) This result was also found by Pond et 
al (1966). The high vaiues cf the water-wave spectra be- 
low 0.05 Hz are not readily explainable. but tay be re- 
lated to a periodic pile-up of water at the leeward end 

of the fetch, where the measurements were taken. Ina 
small lake, such as was used in this study, periodic 

gusts can easily set up a seiche. 

An examination cf the coherence square for the dual- 
probe temperatures shows a high degree of correlation over 
the entire spectrum (Table III). This, together with the 
similarity in the amplitudes of the fiuctuations as seen 


Ne Son View 6LeS then c. 


TQ 


} ato 


in Figure 8, is somewhat surpris 
probe being twice as far froz the water as the lower 

probe. The mean temperatures tor the dual-probe record 
are; 14.26 GC at 2.5 cm and 14,0486 atl Secrpeelnemncam 
air temperature as recorded at approximately 3 ci above 


the water with a mercury therroteter was 14.0°C and the 


2 


Maeer tenperature at the surlLace was 2Q0.33C. Clearly. 
the air near the water was highly unstable, and it ap- 
pears that the greatest air temperature gradient occurs 
Pewew 2.5 cin. 

The results of this study indicate that the data 
collection system used for obtaining terperatures and 
wave records was reasonably accurate and that additional 
work should be done to reduce the level at which the ter- 


PeGaetines Can be taken. 


6. Recommendations for Future Studies 

PUvVEnCEEWworke tiem s tield should inc lude@ageon- 
Pimuous record Of the wind field, additional thermistor 
probes at lower levels (and at least one below the water 
surface), a sinultaneous record of temperature at a fixed 
height above mean water level, and the development of a 
sensor mounting that wili not interfere with the wave 


Field, such as a servo-follower systet.. 
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